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1. INTRODUCTION 

A partition of a natural number n is any integer sequence A = (Ai, A2, . . . ) such 
that Ai > A 2 > • • • > and Ai + A 2 + • • • = n (notation: A h n). In particular, 
Ai = max{Aj G A}. Every partition A h n can be represented geometrically by a 
planar shape called the Young diagram, consisting of n unit cell arranged in consecutive 
columns, containing Ai, A 2 , . . . cells, respectively. 

On the set V n := {A h n} of all partitions of a given n, consider the Plancherel 
measure 

Pn(X):=^, XeV n , (1) 

where d\ is the number of standard tableaux of a given shape A, that is, the total number 
of arrangements of the numbers 1, . . . , n in the cells of the Young diagram A G V n , such 
that the numbers increase in each row (from left to right) and each column (bottom up). 
Note that d\ also equals the dimension of the irreducible (complex) representation of the 
symmetric group & n (i.e., the group of permutations of order n), indexed by the partition 
A (see tZl[TBl EH). According to the RSK 1 correspondence (see |7|), any permutation 
cr G & n is associated with exactly one (ordered) pair of standard tableaux of the same 
shape A G V n . Since there are n\ such permutations, this implies the Burnside identity 
(see ED) 

£4 = < 

thus the measure P n defined in ([TJ determines a probability distribution on V n . 

The Plancherel measure arises naturally in representation-theoretic, combinatorial, 
and probabilistic problems (see J6|). For example, the RSK correspondence implies 



1 Robinson-Schensted-Knuth. 
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that the largest term Ai of the partition A associated with a given permutation a G & n 
equals the length £ n of the longest increasing subsequence contained in a. Therefore, the 
Plancherel distribution of Ai coincides with the distribution of t n in a random (uniformly 
distributed) permutation a E S n (see [2|). 2 

The upper boundary of the Young diagram corresponding to the partition A G V n can 
be viewed as the graph of a stepwise (left-continuous) function X(x),x > 0, defined by 

oo 

X(x) := A l{o}(£) + ^2 = Afa;-], (2) 

i=i 

where 1b{x) is the characteristic function (indicator) of set B and \x] := min{m G 
Z : m > x} is the ceiling integer part of x. Logan and Shepp 1T3l and, independently, 
Vershik and Kerov llT8ll (see also lfl9ll ) have discovered that, as n — > oo, a typical Young 
diagram, suitably scaled, has a "limit shape" represented by the graph of some function 
y = uj{x). This means that for the overwhelming majority of partitions A G V n (with 
respect to the Plancherel measure P n ), the boundary of their scaled Young diagrams is 
contained in an arbitrarily small vicinity of the graph of u(x). 
More specifically, set 

X n (x) := —= X(y/nx), x>0, (3) 
'n 



and consider the function y = cu(x) defined by the parametric equations 

2 

x = - (sm6 - 6cos6), y = x + 2cos6, 0<6<n. (4) 

7T 

The function uj(x) is decreasing on [0,2] and uj(0) = 2, u(2) = 0. Define io(x) as 
zero for all x > 2. Then the random process X n (x) satisfies the following law of large 
numbers IHEBE1: 



We > 0, lim P n < supl X n (x) — u(x) I > e \ 

n-»oo l x >o ' > 



0. 



In particular, for x = it follows that the maximal term in a typical partition asymptoti- 
cally behaves like 2y/n: 

Ai 



> 0, lim Pj — = - 2 > e\ = 0. (5) 

n-»oo I yjn ) 



Remark 1. Due to the invariance of the Plancherel measure under the transposition 
of Young diagrams A <-> A' (when the columns of the diagram A become rows of the 
transposed diagram A' and vice versa), the same law of large numbers holds for X[ = 
#{Aj G A} (i.e., for the number of terms in the random partition A). 



2 The interest in the asymptotic behavior of the random variable l n was stimulated by the Ulam problem 
(see |6 1). The problem was settled by Vershik and Kerov 1 18 1 who showed that l n j ^/n — * 2 in probability 



(cf. formula (|5jl below). 
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A natural question about fluctuations of the random function A n around the limit 
curve uj was posed in [ 1 3 1 (see also |19|), but it remained open for more than 15 years. 
Kerov lfT2l (see also 0) gave a partial answer by establishing the convergence of the 
random process 

A n (x) := y/n (A„(x) - u(x)) (6) 

to a generalized Gaussian process (without any further normalization!). To state this 
result more precisely, it is convenient to pass to the coordinates u = x — y, v = x + y, 
which corresponds to anticlockwise rotation by 45° and dilation by \/2. In the new 
coordinates, the boundary of the scaled Young diagram is determined by the piecewise 
linear (continuous) function \ n (u) and the limit shape is given by 

2 / u 



— (u arcsin — + V A — u 2 ") , \u\ < 2, 

7T V 2 / 

\u\, \u\ > 2. 

Then, according to H21H], the random process 

A n (u) :=v^(An(«)-fi(«)) (V) 

converges in distribution to a generalized Gaussian process A(u), u G [—2, 2], defined 
by the formal random series 

where X 2 , X 3 , . . . are independent random variables with standard normal distribution 
Af(0, 1). Convergence of A n (-) to A(-) is understood in the sense of generalized func- 
tions. It is convenient to choose test functions in the form of the modified Chebyshev 
polynomials of the second kind (see [9|), defined by the formula 

U k (u)\ u=2cose := M{k+ a m , 0€[O,tt]. (9) 

sm 6 

Then one can show (see details in [9]) that for k = 2, 3, ... , 



A n (u) Uk-i(u) du — ► / A(u) Uk-i(u) du 



2 
-2 



~ 9 X 

A(2cosfl) sin(fc^) d6 = ^ 

\Jk 

(here and below, the symbol —>■ denotes convergence in distribution). 
Remark 2. Note that U (u) = 1 (see ©), in which case we have 



A n (u) U (u) du = / A n (u) du = 2 / A n (x)dx 



(10) 
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because the area of the scaled Young diagram equals 1, as well as the area under the 
graph of y = uj{x). This explains why k > 2 in © and (flOb . Also note that fluctuations 
outside [—2, 2] are negligible, since Jj u | >2 <p(u) A n (u) du — > in probability for any test 
function ip with compact support (see ||9|). 

Remark 3. A similar result about convergence to a generalized Gaussian process for 
eigenvalues of random matrices in the Gaussian Unitary Ensemble (GUE) was obtained 
by Johansson iTTOll (see further discussion of these results in [9]). 

However, a "localized" version of the central limit theorem for random partitions 
(i.e., for fluctuations at a given point) has not been known as yet. On the one hand, the 
existence of such a theorem would have seemed quite natural (at least, in the bulk of the 
partition "spectrum" 3 , i.e., for A, G A h n such that i/n ~ x G (0,2)); on the other 
hand, Kerov's result on generalized convergence might cast some doubt on the validity 
of the usual convergence. 

Note that the asymptotic behavior of fluctuations at the upper edge of the limiting 
spectrum (corresponding to x = 0) is different from Gaussian. As was shown in Q for 
Ai and in El [TT1 [HI for any A; with fixed i = 1,2,..., 

Urn pj ^'^f 1 < z) = F t (z), zeR, (11) 

n—*oo I TL I J 

where Fj is the distribution function of the 2-th largest point in the so-called Airy random 
point process, discovered earlier in connection with the limit distribution of the largest 
eigenvalues for random matrices from the GUE (see [ 17 1). In particular, Fi(-) is known 
as the Tracy-Widom distribution function. 

From the point of view of Kerov's limit theorem (see (fluT)). the extreme values 
Ai, A 2 , • • • might present a danger, since according to formula (fTTT) the fluctuations of 
the process A n (x) near x = are large (of order of n 1 ^ 6 ). As this theorem shows, the 
edge of the spectrum in fact does not give any considerable contribution into the inte- 
gral fluctuations. Let us stress, however, that the situation in the bulk of the spectrum 
remained unclear. 



2. MAIN RESULTS 

In our first result, we establish the central limit theorem for the random variable 
A n (x) given by ©. Set 

(.2) 



3 We use the term "spectrum" informally by analogy with the GUE, to refer to the variety of partition's 
terms A; 6 A (cf. the book 1 1 1 where this term is used in a general context of combinatorial structures 
characterized by their components). 
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where 6 X = arccos 



u)(x)— X 
2 



is the value of the parameter 9 in equations © corresponding 



to the coordinates x and y — uj(x). 

Theorem 1. For each < x < 2, the distribution of the random variable Y n (x) with 
respect to the Plancherel measure P n converges, as n — > oo, to the standard normal 
distribution jV(0, 1). 

Remark 4. One can show that Theorem [T] also holds for Y n (x n ) if x n — > x G (0, 2) 

as n — > oo. 

The local structure of correlations of the random process A n (x) is described by the 
following theorem. We write c n x 1 if c„7i, e — > oo, c n n _e — ► for any e > 0, and 

a n x 6„ if a n /b n x 1. 

Theorem 2. Fix xq G (0, 2) and let x±, . . . ,x m G (0, 2) Z?e smc/z ?/za? |xo — %i\ x 
n~ Si l 2 , where < s» < 1 (i = 1, . . . , m). For i — 0, set formally sq = 1. TTzen 
random vector (Y n (x ), . . . , F n (x m )) converges in distribution, as n — > oo, to a Gaus- 
sian vector (Z SQ , . . . , ^ Sm ) w«Y/i zero mean and covariance matrix K with the elements 
K(si, Si) = 1, if (s^ Sj) = min{si, s^} (z 7^ j). 

Note that by Theorem 13 the covariance between Y n (x) and F n (x') asymptotically 
decays as the distance \x — x'\ grows: 



In particular, if \x — x'\ x rr x l 2 (i.e., s = 1), then (Y n (x), Y n (x')) — > Zi), while 
if x' is at a fixed distance from 2 (i.e., s = 0) then Y n (x), Y n (x') are asymptotically 
independent. 

Remark 5. Results similar to Theorems [T] and El were obtained by Gustavsson [ 8 1 
for eigenvalues in the bulk of the spectrum of random matrices in the GUE. 

Let us point out that Theorems [T] and El can be reformulated in coordinates u, v (see 
Sect. 1). To this end, one needs to find the sliding projection (divided by y/2) of the 
deviation A n (u) (see (JT])) onto the line u + v = along the tangent of the graph fi(-) at 
point u. Differentiating equations ©, we get 



x — x x n 



s/2 



=^> lim Cov(Y n (x), Y n (x' 



)) = s. 



n— >oo 



dv x' e + y' e 




du x' e — y' e 



which implies 



A n («) 



26 x 



A n (x n ) (1 + Tl n ) , 



7T 
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where x n — > x and rj n — > (in probability). Hence, setting 

Y n (u):=^M, ueR, (13) 
A/log n 

and using Remark |4[ we obtain the following elegant versions of Theorems [T] and |2j 
where the normalization constant is the same for all points. 

Theorem V . For each —2 < u < 2, the distribution of the random variable Y n (u) 
with respect to the Plancherel measure P n converges, as n — > oo, to the normal distribu- 
tion W(0,1). 

Theorem 2'. Let u , . . . , u m e (—2, 2), \u — Uj| x n^ Si ^ 2 , < < 1, wzY/z r/ie 
conventions as in Theorem^ Then the random vector (Y n (uo), . . . , F„(w m )) converges 
in distribution, as n ^ oo, to a Gaussian vector (Z so , . . . , Z Sm ) with zero mean and the 
same covariance matrix K. 

Remark 6. The covariance function K(s, s') of Theorems 13 and |3 determines a 
Gaussian process Z s on [0, 1] (with zero mean), which can be represented as 

z s = w s + ( s VT^, o<s<i, 

where W s is a standard Wiener process and {( s , < s < 1} is a family of mutually in- 
dependent random variables (also independent of W s ) with normal distribution 7V(0, 1). 
This decomposition shows that the process Z s is highly irregular (e.g., stochastically 
discontinuous everywhere except at s = 1), which is a manifestation of asymptotically 
fast oscillations of the process A n (x) (as well as A n (u)) in the vicinity of each point 
x e (0, 2) (respectively, u G (-2, 2)). 

In conclusion of this section, let us comment on the asymptotics of the random func- 
tion A n (x) at the ends of the limit spectrum, that is, for x = and x = 2. By the 
definition ©, A n (0) = A(0) - ^^(0) = Ai — and according to (fTTT) 

where -Fi(-) is the Tracy-Widom distribution (see Sect. 1). However, the limit distribu- 
tion of A n (2) = A^y 7 ?!) proves to be discrete. 

Theorem 3. Under the Plancherel measure P n , for any z > 

lim P n {A n (2) <z}= F<(0), i = [z] + 1, 

n—*oo 

where Fj(-) is the distribution function of the i-th largest point in the Airy ensemble (see 

dZ2D). 
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Indeed, using the invariance of the measure P n under the transposition A «-> A' (see 
Sect. 1), we have, due to (fTTT) . 

P n {A n (2) <z}= P n {K < 2v^} = Pn ^'J/f 1 < 0} -> *i(0). 
Remark 7. In the "rotated" coordinates it, t>, a similar result holds for both edges: 

lim P n {\ A n (±2) < z} = F f (0), i = [z] + 1. 

n— >oo 

3. POISSONIZATION 

The proof of Theorems [TJand|2]is based on a standard poissonization technique (see, 
e.g., [2 1). Let V = U^ =0 P n be the set of partitions of all natural numbers (as usual, it 
is convenient to include here the case n = 0, where there is just one, "empty" partition 
A(3 h 0). Set |A| := XIa^a ^ an( ^ for t > define the poissonization P l of the measure 



P n as follows: 

2 



F*(A):=e-W^J , AGP. (14) 

Formula (fT4l> defines a probability measure on the set V, since for A G P„ we have 
|A| = n and hence 

°° t n d 2 °° t n 

AeP n=0 ' Ae"Pn ' n=0 

We first prove the "poissonized" versions of Theorems [l] and EJ Let Y t (x) be given 
by formula (fT2"l) with n replaced by t. 

Theorem 4. For each < x < 2, the distribution of the random variable Y t (x) with 
respect to the measure P l converges, as t — > 00, to the standard normal distribution 
Af(0, 1). 

Theorem 5. In the notations of Theorem^ the random vector (Y t (xo), . . . , Y t (x m )) 
converges in distribution (with respect to the measure P l ) to a Gaussian vector (Z so , . . . , 
Z Sm ) with zero mean and the same covariance matrix K. 

In order to derive Theorems [TJ and |2] from Theorems |4] and respectively, one can 
use a standard de-poissonization method. According to formula (fl4l . the expression for 
P t can be viewed as the expectation of the random measure P/v, where N is a Poisson 
random variable with parameter t: 

CO 

P\A) = E [P N (A)\ =e~ t J2u P M)- (15) 

k=0 
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Since N has mean t and standard deviation y/i, equation (TT5b suggests that the asymp- 
totics of the probability P n (A) as n — > oo can be recovered from that of P t (A) as 
£ ~ n — ► oo. More precisely, one can prove that 

P n (A) ~ P*(A), £ ~ 71 -> oo, 

provided that variations of the probability Pk (A) are small in the zone \k—n\ < const-y/n. 
In the context of random partitions, such a de-poissonization lemma was obtained by 
Johansson (see 0). 

4. SKETCH OF THE PROOF OF THEOREM 4 

Note that, in view of ©, the statement of Theorem|4|is equivalent to saying that for 

any zGl 

lim P* [\(Vix) - \Vtx] < 2 V* cos 0s + zy/logi } = $(29^), (16) 

where $(•) is the distribution function of the normal law jV(0, 1). Using the Frobenius 
coordinates Aj — i, set 

P(A):=u^ 1 {A,,-z}, A g P. 

Consider the semi-infinite interval I t := [2\/icos6 x + z\J\og £, oo) and let #I t be the 
number of points \i — i G P(A) contained in J t . Using that the sequence \ — i is 
strictly decreasing and recalling the definition © of the function A(-), it is easy to see 
that relation ([TBI) is reduced to 

lim P<{A G V : #/ t < [v 7 ^]} = $(20^). (17) 

t— >oo 

For = 1, 2, ... , define the A; -point correlation functions by 

Pk{xi, ■ ■ ■ ,x k ) := P*{A G V : xi, . . . ,x k G T>(\)} (xi G Z, Xj 7^ 

The key fact is that the correlation functions pi have a determinantal structure (see 
O): 

Pfc(xi, . . . ,x fc ) = det^Xi,^;*)]^-^, 
with the kernel J of the form 

r ^ j x j v+1 - j x+1 j v ^ x ^ y ^ 

J(x,y;t)=< x-y 

y v£ (j' x J x+ i — J' x+ iJx) , x — y, 

where J m = J m (2y/i) is the Bessel function of integral order m. 
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In this situation, one can apply Soshnikov's theorem [15| (generalizing an earlier 
result by Costin and Lebowitz [5|), stating that the random variable #Jj satisfies the 
central limit theorem: 

v/Var[#/ t J 

provided that Var[#/ t ] — > oo. Thus, in order to derive (fT7t from ([T8l) . it remains to 
obtain the asymptotics of the first two moments of the random variable #Ip The next 
lemma is the main technical (and most difficult) part of the work. 

Lemma 1. Let E f and Var* denote expectation and variance, respectively, under the 
measure P l . Then, as t — ► oo, 

&[#I t ] = VtX-— y/fo£t + 0(l), 
71 

Var*[#J t ] =^(1 + (1)), 

The proof of Lemma 1 is based on a direct asymptotic analysis of the expressions for 
the expectation and variance. In so doing, the calculations are quite laborious and heavily 
use the asymptotics of the Bessel function J m (2y/t) in various regions of variation of the 
parameters. 

Finally, note that the proof of Theorem |5] follows similar ideas using Soshnikov's 
central limit theorem for linear statistics of the form ^ a^I ti ED- 



5. LINK WITH KEROV'S RESULT 

Let us comment on the link between our results and the limit theorem by Kerov lfT2l 
(see Sect. 1). In particular, our goal is to explain heuristically the mechanism of the 
effects that take place for the process A n . 

Note that if \u - u'\ = n~ s l 2 , < s < 1, then 

—2 log | it — it' | 



logn 

That is to say, "time" s indexing the components of the limit Gaussian vector in Theorem 
13, has the meaning of the logarithmic distance between the points u and u', normalized 
by logn. From this point of view, Theorem[2] implies that 



Cov (Y n (u),Y n (u')) w 3 n (u,u') :-- 



mm 



-2 log | it — it' | 
logn 
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In fact, in the course of the proof of Theorems 13 and |2 we obtain that for any e > there 
exist constants Ci, C2 > such that for sufficiently large n, the following estimate holds 
uniformly in u, v! E (—2, 2): 

Cov(A„ (lI ),A„ M )<( "?'° g !""l! e ' < 19 > 

I C 2 , \u- u'\ > e. 

Consider now the integral of A n (u) with respect to a test function (p: 

A„M := J A n (u) <f(u) du. 

Using (IT9l) we have 

Var (A n [cp]) = / / (p(u) (p(u') Cov(A n (u), A n (u')) du du' 
J -2 2 

< — Ci // </?(it) y(u') log I it — it' I <itt dw' 

J J\u-u'\<e 



+ C 2 / / ip(u) (p(u r ) dudu' < 00, 

J J\u-u'\>e 

since the function log |it — it'| is integrable at zero. Therefore, A„ [<^] is bounded in 
distribution as n — > 00, which helps understand why Kerov's result holds without any 
normalization (see Sect. 1). 

Remark 8. We believe that by sharpening the asymptotic estimates (fT9l) . it may be 
feasible to compute the limit of the variance Var(A n [<£>]) and thus recover the Kerov 
theorem directly from the analysis of the correlation structure. We will address this issue 
elsewhere. 

Conversely, the limiting process A (it) defined in ® can be used to get the informa- 
tion contained in Theorems [T] and|2 (at least heuristically). To this end, observe that 
since the number of terms in a typical partition is close to 2y/n (see Remark [T]), it is 
reasonable to think that the number of "degrees of freedom" of a random partition A h n 
is of order of m x ^fn, and hence the random variable A n (it), u E (—2, 2), may be 
represented by the partial sum of the series ® 

n / m 2-^\X k sm(k9) . . 

Note that for any u = 2 cos 9, u' = 2 cos 9' 

m 



7T 2 ^— ' k 

k=2 

2 j\ cos{k{9 - 9')) 2 y> cos{k{9 + 9')) 

fc=2 k=2 
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The second sum in (l20b converges for all 9, 9' E (0, 7r). For 9 = 9', from d20b we get 

2 Al 21ogm logra 
Var[S m ( M )]~-^-~ — — 

k=2 

and it follows that 



f. — ►A (0,1) (n->oo), 

Vlogra 

which is in agreement with Theorem[l]. Moreover, if u' — u x n~ s l 2 (and hence 9' — 9 
77,- s / 2 ) then the first sum in (|2~0"1) is approximated by the integral 

cos(x(e'-e)) , r m|9 '- e| cosy , r 1 , 



Hence 



x J2\o'-e\ V J2\e'-e\ V 

s 
2 



— log 1 6' — 9 1 ~ - log n. 



Covf.S'm^), SV^m')) ~ ^ (n -> oo), 
as predicted by Theorem EI- 

Remark 9. As already mentioned (see Remarks |3] and |5] and also a comment after 
formula <HU)), there is similarity between the asymptotic properties of the spectra of 
random partitions and random matrices from the GUE. Our discussion suggests that the 
relationship between the generalized type convergence ifTOll and the localized central 
limit theorem [8] in the GUE can also be explained using the correlation structure of the 
spectrum. One can expect that similar relationship may be in place for other classes of 
random matrices and for more general determinantal random ensembles, but this issue 
needs to be studied further. 

Remark 10. Let Y(u), u e [—2, 2], be a random process with independent values, 
such that Y{u) has a standard normal distribution for each u E (—2, 2), and ^(±2) = 
(a.s.). Our results (see TheoremEI and Remark|7] after Theorem^ imply that the random 
process Y n (-) (see (fT3l) N ) converges to Y(-) in the sense of finite dimensional distributions. 
A natural question may arise as to whether this can be extended to weak convergence. 
However, it is easy to see that the answer is negative, at least under the natural choice of 
the space of continuous functions C[— 2, 2], because the necessary condition of tightness 
breaks down (see Q). Indeed, for any 5 > 0, e > and all u, u' E (—2, 2) such that 
| it — u'\ < 5, we have 

lim P n {\Y n (u) - Y n {u')\ >e}= P{\Y(u) - Y{u')\ > s) > 0. 

n— >oo 

Analogous remark applies to the process Y n (x), x E [0,2], considered in the space 
D[0, 2] of right continuous functions with left limits. 
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